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SECTION  I 
INTRODUCTION 


BACKGROUND 

The  Naval  Surface  Weapons  Center  (NSWC)  is  attempting  to  obtain  a 
practical  computer  solution  to  the  air  and  water  shock  wave  propagation  re- 
sulting from  a high  explosive  (HE)  detonation  above  the  surface  of  the  water. 
Earlier  hydrocode  calculations  of  the  HE  problem  were  done  using  an  Eulerian 
wave  propagation  code  and  an  Eulerian  code  coupled  with  an  acoustic  wave 
propagation  code.  The  results  from  one  such  calculation  reported  In  reference 
1,  Indicate  that  the  Eulerlan/acoustic  shock  wave  calculation  generated 
reasonable  results.  However,  the  resolution  of  the  shock  wave  Is  not  as  good 
as  could  be  expected  from  Lagranglan  shock  wave  calculations  where  the  grid 
moves  with  the  material  and  material  Interfaces  are  maintained.  Therefore, 
the  NSWC  Is  attempting  to  obtain  a practical  solution  using  a Langranglan  wave 
propagation  code. 

The  two  codes  proposed  for  feasible  use  were  TOODY  II  (ref.  2)  and  STEALTH 
(ref.  3).  After  Investigation  Into  both  codes  it  was  decided  that  STEALTH 
would  require  modification  beyond  the  scope  of  this  task  to  handle  the  size 
and  difficulty  of  the  problem  proposed.  Also,  the  automatic  rezoning  capability 
was  not  general  enough  to  handle  the  severity  of  deformations  expected  In  this 
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Cylinder  Above  a Water  Surface:  an  Interim  Report , N5WC/W0L/TR-77-34, 

Naval  Surface  Weapons  Center,  White  Oak  Laboratory,  Silver  Spring,  MD, 

January  1978. 

2.  Bertholf,  L.D.,  Benzley,  S.E.,  TOODY  II,  A Computer  Program  for  Two-Dimensional 
Ware  Propagation,  SC-RR-68-41,  Sandia  Laboratories,  Albuquerque,  New  Mexico, 
November  1968. 

3.  Ronald  Hofmann,  et  al , STEALTH , A Lagrange  Explicit  Finite-Difference  Code 
for  Solids,  Structural,  and  Thermohydraulic  Analysis,  EPRI-NP-160, 
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problem.  The  author's  extensive  experience  with  TOODY  II  In  calculating  wave 
propagation  In  geologic  materials  and  the  generality  of  the  rezoning  code 
TOOREZ  (ref.  4)  were  also  factors  In  the  selection  of  these  codes  for  the 
calculation. 

OBJECTIVE 

The  purpose  of  this  effort  was  to  verify  that  a Lagranglan  wave  propaga- 
tion code,  namely  TOODY  II,  could  be  used  to  obtain  a solution  to  the  HE 
detonation  over  water  problem,  and  that  a higher  resolution  of  the  water  shock 
wave  would  be  obtained  by  using  the  Lagranglan  frame  of  reference. 

SCOPE 

This  effort  concentrated  on  obtaining  a solution  to  the  HE  detonation 
over  water  problem  In  the  form  of  pressure  time  histories  at  the  air/water 
Interface  and  at  various  depths  along  the  axis  at  symmetry. 


4.  Thorne,  B.J.,  and  Holdrldge,  The  TOOREZ  Lagrangian  Rezoning  Code , SLA- 
73-1057,  Sandla  Laboratories,  Albuquerque,  NM,  April  1974. 
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SECTION  II 
PROBLEM  DESCRIPTION 


INITIAL  CONDITIONS 

The  general  problem  for  which  a practical  solution  is  desired  consists 
of  a 454  gm  (1  lb)  cylinder  of  pentolite,  17  cm  long  and  4,6  cm  tn  diameter 
Ignited  in  air.  The  center  of  the  cylinder  is  located  from  15.24  to  30.48  cm 
(.5  to  1.  ft)  above  the  water  surface.  For  the  specific  problem  discussed  in 
this  report,  the  center  is  located  at  30.48  cm  (1  ft)  above  the  water.  This 
problem  was  selected  because  the  expected  pressure  in  the  air  shock  upon 
impact  on  the  water  would  be  less  than  that  for  the  15.24  cm  height-of-burst. 
This  would  simplify  the  calculation  during  the  air  shock  transmission  into  and 
the  reflection  off  of  the  water.  The  problem  is  illustrated  in  figure  1. 


The  problem  has  an  axis-of-symmetry  down  the  center  of  the  penolite  cy- 
linder and  a plane-of-symmetry  at  the  mid-height  of  the  cylinder  if  the  water 
Is  Ignored.  Therefore,  only  a quarter  of  the  problem  had  to  be  modelled  in 
the  calculation,  also  shown  in  figure  1. 


The  equation-of-state  (EOS)  used  for  the  detonation  gas  products  was 
the  JWL  equation  for  pentqlite,  The  equation  for  pressure  is  as  follows: 


P = A 1.0  - 


R*Po 


_Ri  Pjl 

» p + B 1 1.0 


WiL_ 

R2Po 


Rz.fLSL 


+ VJps 


where. 


p = density 
e = energy  density 


{ 
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A = 5.3177 
B = 0.08933 
Ri  = 4.6 
R2  = 1.05 
W = 0.33 

And  the  initial  conditions  are, 

D0  = .736  cm/ysec  (detonation  velocity) 

b5  = 2.5,  constant  determining  width  of  detonation  wave 

p0  = 1-65  gm/cm3 

e0  = .080  Mbar-cm/gm  (reference  energy  density) 

The  air  was  modelled  as  an  ideal  gamma-law  gas  where  pressure  is  calculated 
as  follows: 

P - (r  - 1)  pe 

where , 

T = 1.4 

and  the  initial  conditions  are, 

Po  = .001225  gm/cm3 
C0  = .0340294  cm/ysec 

A hydrodynamic  EOS  was  used  to  model  the  water  where  pressure  is  calculated 
as  follows: 

P = f i (p)  + f2(p)e 

where , 

f,  - p„(  i - Jt] 

f2  = rP 


5 


PH  = Kon(  1 + + k2H  + K3n  ) 


r = r0  (l-n) 

E = 1 


Constant  Tp 


y •£--  1 

Po 

K0  a PoCo  = .0219504 
Ki  = 3.7815448 
K2  = 2.8199746 
K3  = 115.79338 

and  the  initial  conditions  are, 

Po  = .9982  gm/cm3 
C0  = .14785468  cm/ysec 
Pmin  = -l.xlO-6  Mbar  (1  atm  tension  cutoff) 
v = .5  (Poisson's  ratio,  hydrodynamic) 

For  details  of  the  equations  of  state  used  in  TOODY  II  see  ref.  3. 


The  initial  finite  difference  grid  for  the  calculation  was  setup  for 
efficient  core  storage  utilization  and  for  ease  of  future  rezoning  by 
defining  the  overall  grid  size,  locating  the  air  water  interface  and  fine 
zoning  the  active  region  of  the  problem  only.  The  HE  and  air  region  nearby 
were  zoned  with  approximately  .25  by  .25  cm  zones.  Outside  this  region 
zones  extended  to  edge  of  the  problem  at  a range  of  200  cm,  to  the  air/water 
interface  at  30.48  cm  and  from  the  air/water  interface  to  a depth  of  70  cm. 

This  minimizes  the  number  of  calculations  that  have  to  be  done  during  any 
given  time/cycle  by  minimizing  the  number  of  grid  points.  The  grid  at  the 
beginning  of  the  calculation  is  shown  in  figure  2 and  details  of  the  HE 
zoning  are  shown  in  figure  3.  The  grid  plots  in  the  remainder  of  the  report  show 


Figure  2.  Initial  Grid  for  the  Lagrangi 


only  the  grid  in  the  region  of  the  shock  wave  (the  region  where  pressure  has 
exceeded  1 bar).  It  should  be  remembered  that  the  complete  grid  extends  to  a 
range  of  200  cm  and  a depth  of  70  cm. 

As  the  calculation  progressed  zone  sizes  were  selected  to  maintain  a 
reasonably  uniform  grid.  The  air  and  HE  zones  varied  from  1.5  by  2.5  cm  to 
l.by  .5  cm.  When  the  air  shock  reached  the  surface  of  the  water  the  water 
was  zoned  as  shown  in  figure  4. 

REZONING  PROCEDURE 

In  a Lagrangian  frame  work  the  reference  grid  moves  with  the  material  so 
that  in  the  region  of  shock  fronts  a grid  can  become  severely  distorted  par- 
ticularly for  a highly  compressible  material  like  air.  As  the  shock  wave 
propagates  through  the  air  it  crushes  up  air  zones  resulting  in  a very  small 
timestep  and  possibly  unstable  grid  configurations.  Also  as  the  HE  gas  pro- 
ducts approach  the  water  surface  air  zones  are  compressed  and  the  air  is 
forced  to  flow  out  from  between  the  water  and  explosive  products.  For  these 
reasons  rezoning  is  required  to  redistribute  the  grid  points  in  the  region  of 
shocks,  to  remove  air  cells  from  between  the  HE  gas  products  and  the  water  and 
to  place  a frictionless  sliding  interface  at  the  air/water  interface. 

Three  criteria  were  established  for  determining  when  the  grid  needed  to 
be  rezoned.  They  were  (1)  when  the  timestep  dropped  by  more  than  a factor 
of  2 in  any  one  cycle  indicating  an  unstable  zone,  (2)  when  the  timestep  drop- 
ped below  some  selected  minimum  timestep  or  (3)  when  the  activity  was  within 
one  zone  of  the  edge  of  the  well  defined  or  finely  zoned  region  of  the  prob- 
lem. 


The  general  procedure  for  the  T00DY/T0RREZ  calculations  consisted  of  a 
T00DY  execution,  a restart  file  generation  of  the  deformed  grid,  a rezoning 
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Air  and  HE  grid  determined 


Figure  4.  Initial  Zoning  Definition  for  the  Water  as  Generated  by  Rezoning  at 
59.6,  62.1  and  66.9  Microseconds 


of  the  deformed  grid,  and  restarting  of  the  TOODY  calculation  from  the 
rezoned  restart  file.  For  this  calculation  this  procedure  was  repeated  23 
times  to  a problem  time  of  113  psec. 

Figures  5,  6 and  7 show  deformed  and  rezoned  grids  at  three  significant 
times.  Shown  in  figure  5A  is  a grid  plot  at  10.7  psec.  It  shows  that  about 
75%  of  the  pentol ite  has. detonated  and  that  the  shock  wave  is  just  starting 
to  pull  away  from  the  HE  products  at  the  upper  symmetry  boundary.  Figure  5B 
is  the  rezoned  grid  at  the  same  time  showing  the  redistributed  grid  points  in 
the  shock  wave  and  that  the  HE/air  interface  has  been  maintained. 


Figure  6A  shows  the  grid  at  59.6  psec.  The  entire  pentolite  cylinder 
has  been  detonated  and  the  almost  spherical  shock  front  is  evident  by  the 
region  of  very  dense  zoning.  The  crushing  up  of  the  air  cells  between  the  HE 
products  and  water  is  also  quite  evident.  Figure  6B  shows  the  rezoned 
grid  at  the  same  time.  The  figure  shows  that  the  grid  points  in  the  shock 
front  have  been  redistributed,  air  cells  haye  been  removed  from  between 
HE  products  and  the  water  and  very  fine  zoning  of  the  water  has  been  added. 

Shown  in  figure  7 is  the  deformed  grid  at  the  termination  of  the  problem 
at  113  psec.  It  shows  the  shock,  wave,  to  be  ahout  5 cm  out  in  front  of  the 
HE  products.  It  has  also  propagated  about  24  cm  along  the  surface  of  the  water. 
Also  evident  in  the  figure  is  the  reflected  shock  wave  located  5 cm  above  the 
water.  The  calculation  was  terminated  at  this  point  and  therefore,  this 
grid  was  not  rezoned. 


11 


d at  59.6  Microseconds  (Before  Rezoning) 


croseconds  (After  Rezoning) 


SECTION  III 


RESULTS  OF  THE  CALCULATION 


The  calculation  was  terminated  at  a problem  time  of  113  usee.  At  this 
time  the  shock  wave  had  propagated  about  7.5  cm  Into  the  water.  The  results 
of  the  calculation.  In  the  form  of  pressure  and  impulse  time  histories  at 
a range  0.0  cm  (along  the  symmetry  axis)  and  depths  of  0.0  and  4.2  cm,  are 
shown  In  figures  8,  9 and  10.  There  are  two  sets  of  data  at  the  air/water 
Interface.  One  set  (designated  -0.  cm)  is  the  pressure  and  impulse  in  the 
air  adjacent  to  the  water  surface.  The  second  set  (designated  +0.  cm)  is  the 
pressure  and  Impulse  In  the  water  cell  just  below  the  Interface. 


Other  Important  numerical  results  from  the  computation  are  the  "shock 
factors"  which  are  semi -empirical  quantities  that  correlate  well  with  target 
damage.  It  Is  defined  as  follows: 

Q » W/R 

where, 

W - charge  weight  in  lbs  of  TNT 
R ■ range  In  feet  from  the  explosive 

It  has  been  determined  that  an  "equivalent  shock  factor"  resulting 
from  an  alrblast  can  be  calculated  from  the  analytical  expression  that 
follows: 


Qe 


Ml 

A P‘ 

— ITPT5 — 


.02416^  1 


iMfli 


Pressure  and  Impulse  at  -0.0  CM  (Air  Pressure) 


Ill 


Figure  10.  Pressure  and  Impulse  at  4.2  CM 


t 


P = peak  pressure  in,  psi 

I = impulse  in  psi-msec 

E = energy  flux  density  in  in-lbs/in2 


* / P(t)2dt 
PoC0 

For  the  pressure  and  impulse  time  histories  in  figures  8,  9 and  10  the 
equivalent  shock  factors,  Qe,  are  shown  in  Table  1 along  with  pertinent  peak 
parameters  used  in  Qe  calculations. 


* 

) 


DEPTH, 

CM 


TABLE  1 

EQUIVALENT  SHOCK  FACTORS  AND  PEAK  PARAMETERS 


PEAK 

PRESSURE, 

BARS 


PEAK 

PRESSURE, 

PSI 


7194. 


IMPULSE, 

ENERGY 

SHOCK 

FLUX, 

FACTOR 

PSI-MSEC 

IN-LBS/IN2 

Qe  ** 

138. 

185. 

.28 

128. 

160. 

.25 

90* 

60* 

.15* 

* based  on  linear  extrapolation  of  pressure 
' to  zero  in  50  ysec 

k* 

Qe  - equivalent  shock  factor  for  454  gm  pentolite  cylinder 
at  30.48  cm  HOB 


SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  calculations  presented  in  this  report  were  not  of  long  enough 
duration  for  direct  comparison  with  or  even  extrapolation  to  experimental 
results.  However,  an  assessment  of  the  accuracy  can  be  made  by  comparisons 
with  previous  Eulerian  and  Eulerian/acoustic  code  calculations  of  the  same 
problem.  Preliminary  comparisons  indicate  that  at  a minimum  better  resolution 
of  the  shock  wave  was  attained. 

Several  problems  encountered  during  the  calculation  could  be  better 
handled  or  avoided  in  future  attempts  at  obtaining  a solution  via  Lagrangian 
wave  propagation  codes.  The  major  recommendation  would  be  to  obtain  a new 
code  or  modify  an  existing  one  incorporating  a problem  oriented  automatic 
rezoning  capability.  Version  of  TOODY  and  STEALTH  have  automatic  rezoning 
features  that  might  be  modified  to  handle  this  type  of  problem.  However, 
a major  effort  is  probably  required. 

If  the  procedure  used  for  this  calculation  is  repeated  or  continued, 
special  caution  should  be  used  when  the  shock  wave  comes  in  contact  with  the 
water.  Particularly  away  from  the  symmetry  axis  where  the  air  shock  has 
other  then  normal  incidence  and  during  the  relief  portion  of  the  loading.  Un- 
reasonably large  tension  can  be  generated  due  to  the  zoning  configurations 
on  either  side  of  the  air/water  interface.  Investigation  is  required  into  what 
combinations  of  zone  sizes  will  facilitate  accurate  transmission  and  reflection 
of  the  air  shock  wave.  Zone  sizes,  aspect  ratios  and  artificial  viscosity 
coefficients  have  been  shown  to  have  very  significant  effects  on  Lagrangian  wave 
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propagation  calculations  as  discussed  in  reference  6. 


This  report  verifies  that  a Langragian  wave  propagation  code  can  be  used 
to  obtain  a reasonably  accurate  solution  to  the  high  explosive  detonation  over 
water  problem.  However,  the  number  of  rezonings  required  (23  between  0 and 
113  ysec)  and  the  awkwardness  of  multiple  rezonings  and  restartings  leave  a 
judgement  on  the  practicality  and  accuracy  of  the  procedure  outlined  herein 
open  to  debate. 
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5.  Higgins,  C.J.,  Rudeen,  D.K.,  Effects  of  Zone  Size , Aspect  Ratio , and 
Artificial  Viscosity  on  the  Accuracy  of  Ground  Motion  Calculations  with 
Finite-Difference  Codes,  AFWL-TR-77-108,  Air  Force  Weapons  Laboratory, 
Kirtland  AFB,  New  Mexico,  September  1977. 


t 


24 


SECTION  V 


REFERENCES 


1.  Lutsky,  M. , Eulerian  Hydrocode  Calculations  of  the  Detonation  of  an 
Explosive  Cylinder  Above  a Water  Surface:  an  Interim  Report , NSWC/WOL/ 
TR- 77-34,  Naval  Surface  Weapons  Center,  White  Oak  Laboratories,  Silver 
Spring,  MD,  January  1978. 

2.  Bertholf,  L.D.,  Benzley,  S.E.,  TOODY  II,  A Computer  Program  for  Two- 
Dimensional  Wave  Propagation , SC-RR-68-41,  Sandia  Laboratories, 
Albuquerque,  New  Mexico,  November  1968. 

3.  Ronald  Hofmann,  et  al , STEALTH,  A Lagrange  Explicit  Finite-Difference 
Code  for  Solids,  Structural,  and  Thermohydraulic  Analysis,  EPRI-NP- 
160,  Electrical  Power  Research  Institute,  Palo  Alto,  CA,  August  1976. 

4.  Thorne,  B.J.,  and  Holdridge,  The  TOOREZ  Lagrangian  Rezoning  Code, 
SLA-73-1057,  Sandia  Laboratories,  Albuquerque,  New  Mexico,  April  1974. 

5.  Higgins,  C.J.,  Rudeen,  D.K.,  Effects  of  Zone  Size,  Aspect  Ratio,  and 
Artificial  Viscosity  on  the  Accuracy  of  Ground  Motion  Calculations  with 
Finite-Difference  Codes,  AFWL-TR-77-108,  Air  Force  Weapons  Laboratory, 
Kirtland  AFB.  New  Mexico,  September  1977. 


25 


